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Abstract. We provide a unified description of the hemispherical asymmetry in the cosmic 
microwave background generated by the mechanism proposed by Erickcek, Kamionkowski, 
and Carroll, using a 5J\f formalism that consistently accounts for the asymmetry-generating 
mode throughout. We derive a general form for the power spectrum which explicitly exhibits 
the broken translational invariance. This can be directly compared to cosmic microwave 
background observables, including the observed quadrupole and /nl values, automatically 
incorporating the Grishchuk-Zel’dovich effect. Our calculation unifies and extends previ¬ 
ous calculations in the literature, in particular giving the full dependence of observables on 
the phase of our location in the super-horizon mode that generates the asymmetry. We 
demonstrate how the apparently different results obtained by previous authors arise as dif¬ 
ferent limiting cases. We confirm the existence of non-linear contributions to the microwave 
background quadrupole from the super-horizon mode identified by Erickcek et al. and fur¬ 
ther explored by Kanno et al., and show that those contributions are always significant in 
parameter regimes capable of explaining the observed asymmetry. We indicate example pa¬ 
rameter values capable of explaining the observed power asymmetry without violating other 
observational bounds. 
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1 Introduction 

There is substantial variety of work in the literature proposing mechanisms behind the hemi¬ 
spherical anomaly in cosmic microwave background (CMB) temperature first identified in 
WMAP data [1], and later corroborated by the Planck satellite collaboration [2] and others 
[3-5]. We take the observed value of the power asymmetry to be A = 0.06 [6], where 2 A 
is the amplitude of the power asymmetry in the direction where it is maximal. The most 
popular scenario is that proposed by Erickcek et al. [7, 8], which we refer to as the EKC 
mechanism. This proposes that the power spectrum is modulated due to the presence of a 
high-amplitude long-wavelength perturbation, the wavelength being substantially larger than 
the present observable Universe. This scenario is not viable for perturbations generated by 
single-field inflaton, but can be if an auxiliary field such as a curvaton is responsible. 

More generally, cosmological models that give rise to high-amplitude field fluctuations 
on long wavelengths can imprint distinct signatures on the small-scale density perturbations, 
by breaking the translational invariance of the background in our observable patch. This in 
turn suggests the possibility of probing scales beyond our observable Universe by searching 
for signs of broken translational invariance in the sky. A particular example scenario with 
such large fluctuations on long wavelengths is proposed in [9] , which identifies the large-scale 
mode with a super-curvature mode excited by vacuum tunnelling, providing for the first time 
a physical framework in which complete calculations can be undertaken to test against the 
observed properties of perturbations including their asymmetry. 

However a unifying formalism for describing the hemispherical anomaly is still missing, 
one that holds across different treatments in the literature, and describes the full set of effects 
that it imparts on CMB observables. Most of the results assume specific regimes where a 
given set of assumptions holds. For this reason it is not always clear how results from different 
models can be compared, without repeating the calculation anew under the new assumptions, 
with the consequence that results in the literature don’t always readily carry across different 
works. 
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In this study we start with a field perturbation which is a combination of the usual 
continuum inflationary perturbation spectrum together with a super-horizon mode; the latter 
is typically inserted by hand though it may have a physical origin as a super-curvature mode 
also generated by inflation. We proceed by treating the resulting curvature perturbation 
‘blindly’ using the SAf formalism, i.e. we compute the power spectrum with a minimal set 
of assumptions and retain all contributions of the super-horizon mode including non-linear 
terms and interactions between the short- and long-wavelength modes. 

Other studies in the literature also considered a global perturbation composed of the two 
perturbation components [7, 8, 10-16], some of which invoked the 5J\f formalism, however 
have not carried out an exhaustive computation keeping all linear and non-linear terms. 
This has led to some confusion in the literature regarding what observational signals are 
actually induced by super-horizon Held fluctuations. We clear up the disagreements among 
the literature and unify the findings of previous studies by consistently incorporating all 
fluctuations in the 5AT calculations. In particular, we confirm the existence of non-linear 
contributions to the quadrupole first identified in Ref. [8] and further analyzed by Kanno et 
al. [14], then extend those studies to include higher-order terms. The authors of the EKC 
mechanism did not include non-linear corrections to the power asymmetry or non-Gaussianity 
parameter /nl in Ref. [7], but included these for their quadrupole expression in Ref. [8]. In 
the regime of small non-linearities where their calculations are fully valid, the amplitude of 
the generated asymmetry is too small to fit the observed value, and so cannot be responsible 
for the hemispherical anomaly, whereas our results are valid in the regime where sufficient 
asymmetry can be obtained. 

Moreover, our unified treatment yields the full shape of the power spectrum, not dis¬ 
tinguishing between super-horizon and continuum mode contributions to multipoles. This 
allows us to compare our values directly to CMB data. In particular we can compare the 
derived quadrupole and octupole with their observational bounds, avoiding separate eval¬ 
uations of the Grishchuk-Zel’dovich (GZ) [17] and EKC effects on the quadrupole as in 
Ref. [11, 14]. There has been discussion in the literature as to whether these contributions 
are to be regarded as correlated or independent, particularly because the GZ effect implies 
a stochastic treatment of the super-horizon mode while the EKC mechanism presupposes an 
asymmetry with a well-defined direction and phase. For this reason, in our formalism we 
treat the large-scale mode non-stochastically. The observed direction of the asymmetry is 
generated by a single realization of the large-scale mode and associated to a specific phase, 
and we can then no longer consider its phase and direction to be random. In particular we 
give the full dependence of the different observables on the phase which the observer occupies 
on the super-horizon mode. 1 

Lastly our expressions for observables in the asymmetric sky significantly simplify those 
in the literature by identifying the correct background field values for evaluation. We obtain 
observables like the power spectrum and the non-linearity parameter, /nl> by keeping their 
usual expressions when no anomaly is present, but evaluating them at a shifted value of the 
field. We find that shifting the homogeneous background value of the field to the value where 
we are in the large-scale perturbation correctly accounts for the effects that the asymmetric 
mode imparts on the sky. With this simplification we obtain the correct expressions for 
observables, including non-linear effects, while using standard expressions valid when no 
asymmetry is present. 

1 Previous works such as [7, 8] also presented the phase dependence for some quantities. 
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In Sections 2 to 5 we use the formalism described above, treating the unified field 
perturbation as the source of the overall observables and keeping all terms up to second order 
in the field derivatives. In the regime where second-order terms are dominant, considered by 
Kanno et al. [14], we arrive at similar expressions to theirs. However, having found physical 
regimes where non-linear terms dominate we then have reason to investigate even higher- 
order terms as well, and we proceed to do that in Section 6. There we arrive at the most 
general expressions for the effects the asymmetry imparts on the power spectrum, bispectrum 
and trispectrum, with full inclusion of non-linear terms and interactions between short- and 
long-wavelength modes. These expressions are new to the literature. 


2 Power spectrum with broken translational invariance 


In this article we use the 5J\T formalism [18-21] to describe the hemispherical asymmetry. 
The curvature perturbation £ is generated by fluctuations in a field a, typically a curvaton 
field in the context of the EKC mechanism [7, 8]. We are going to split the fluctuations of the 
field a into two separate contributions 5a and A a; later on we will identify 5a as the ‘normal’ 
continuum spectrum of fluctuations sourced by the inflationary background, and use Act to 
represent the super-horizon field fluctuation invoked by the EKC mechanism. The curvature 
perturbation £ is determined by the variation in the number of e-foldings at different spatial 
locations and can be expanded in terms of the field perturbation as 

£(®) = A f a ( 5a(x ) + Aa(x)) + ^A/^ (5a{x) + Acr(a;)) 2 4-. (2.1) 

Here A f is the number of e-folds counted forward in time from an initial flat hypersurface to 
a final uniform-density hypersurface, and N a , Afua, • • • denote derivatives of M in terms of 
the value of a on the initial hypersurface. The arguments of Af a , Af„ a , ■ ■ • are taken as the 
background field value <7b g of the entire universe in the expression (2.1). Hereafter, unless 
explicitly stated otherwise, the derivatives of M are evaluated at the background value Cbg. 
We remark that on some very large scale the Universe is assumed to be homogeneous and 
isotropic, so that the separate universe picture is valid. 

Going to the Fourier space via 

C(*) = (2^)3 / Cfc ’ (2 ' 2) 


one obtains 


Cfc = ACx (5a k + Aa k ) 4- j d 3 p ( 5a p + Acj p ) (5a k ~ p + Aa k _ p ) 4-• (2.3) 


To compute the power spectrum, we assume 5a and A a to be uncorrelated to each other, 
and also that the small-scale field perturbations 5a are well approximated as Gaussian, so 
that 

{5a k ) = 0, (2.4) 

(5a kl 5a k2 ) = (2 tt) 3 P s<J {ki)5(k 1 + fe 2 ), (2.5) 

(5a kl 5a k2 5a k 3 ) = 0, (2.6) 

where k = |fc|. Even though we suppose the bispectrum of 5a to vanish, non-Gaussianities 
in £ will be generated via the non-linear relation between £ and 5a. We will see that such 
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local-type non-Gaussianities are directly related to the hemispherical asymmetry. One can 
also study cases with other types of non-Gaussianities arising directly from the bispectrum 
of the field fluctuations themselves, by modifying the assumption ( 2 . 6 ). 

Furthermore, we suppose Aa k to vanish at the wavenumbers on which the power spec¬ 
trum will be measured. In other words, Act k is assumed to be nonzero only at wavenumbers 
much smaller than the usually considered wavemodes. Then, even if Aa k is stochastically 
generated, we can only observe one realization in our current Hubble patch and thus we will 
treat A 07 , as a non-stochastic quantity in the following calculations . 2 Then we obtain 

(CfciCfc 2 ) =( 27 r) 3 A /' 2 Psa{k 1) 6(ki + k 2 ) + Af a J\T aa Aa kl+k2 (Psa(h) + Psa(k 2 )) 

a r 2 r , 

+ 777 73 / d pAa p Aak 1 +k 2 -pPsc 7 (|^i p\) 

w f 3 3 (2 - 7) 

+ J d pd' q (Aa p Acr q A(T kl ~ p Acr k2 _ q + ( 5a p 5cr q 5a kl _ p 5a k2 - q )) 

+ (terms with J\f aacr , N aaa<y , ■■■)■ 

Here terms that involve 5{k\) and 5(k 2 ) were dropped, as we do not consider power spectra 
at zero wavenumber. 

For the super-horizon fluctuation Act, let us assume that it only has modes at a single 
wavenumber Aki,, as is often assumed and for instance is the case in the open Universe 
scenario proposed in Ref. [9]. Then we can write 

Aa k = ^Y~Aa{5(k - fc L ) + 5(k + k L )} , (2.8) 

where Aa is real and denotes the fluctuation amplitude. Two delta functions are introduced 
so that the fluctuation in real space is real-valued, 

A a(x) = A a cos ( k l • x ). (2.9) 

Here it should be noted that A a has no spatial dependence, but can depend on time. Then, 
taking into account that the wavenumbers under consideration for the power spectrum are 
much larger than &l, we arrive at 

( 2 7 r ) 3 (Cfc i Cfc 2 ) = Af^Psaik^Siki + k 2 ) 

+ Af a KaP S a(kl) Aff {<*(*! + k 2 - k L ) + (J(fc! + k 2 + k L )} ( 2 . 10 ) 

+ — Afp a Pscr(ki) Acr~ {25(fci + k 2 ) + 5(ki + k 2 — 2&l) + ^(^i + k 2 + 2k l)} . 

Here we have assumed that P§ a is nearly scale-invariant, i.e., 

Psa (\ki ± k L \) = Ps a (fci). ( 2 -ll) 

2 Generally speaking, the large-scale fluctuation Act can induce inhomogeneities in the power spectrum 
of the small-scale fluctuations 5a through the field’s self-interactions, giving rise to translational-invariance 
breaking in the right-hand side of Eq. (2.5). Furthermore, the bispectrum can evolve in time and thus its value 
does not necessarily stay zero as in Eq. (2.6). For instance, in cases where /nl is strongly scale-dependent, 
the time evolution of the correlation functions of 5a may become significant. However, for simplicity, in this 
article we neglect such effects and assume the power spectrum and bispectrum to take the forms (2.5) and 

(2.6) at all times. 
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We have also neglected the (<5<r 4 ) term in Eq. (2.7) as 5a is small and its effect will be 
negligible unless /nl is huge. Furthermore, terms with Ak a aa, A! aa aa , ■■■ were dropped. 
Such higher-order derivatives of AT are neglected hereafter, until Section 6 where we provide 
an alternative and more general calculation. 

The second and third lines in Eq. (2.10) show that the power spectrum of ( does not 
obey translational invariance, because of the presence of the super-horizon fluctuation A a. 
In the limit ki, —> 0, the power spectrum reduces to 

(CfciOfc 2 > = (27T) 3 P 5a (k 1 )5(k 1 + k 2 ) (Af a +Af aa Aa) 2 (2.12) 

= (2ir) 3 P$ a (ki)5(ki + fc 2 ) (A/^^-) 2 , (2.13) 

where in the second line we have absorbed A a into the argument of AT a . Thus it is clear 
that in this limit the super-horizon fluctuation is simply a homogeneous shift Act of the 
background field value. The Taylor series in Eq. (2.12) is taken only up to quadratic order 

in Act, because we have been neglecting third and higher derivatives of AT in the calculations. 

The asymmetric nature of the power spectrum is clearer when moving to real-space: 


K0>0C(»)> = (Tp f 

3 1 . [ Ipke' k Wy)p 6(7 (fc) {ACr + M aa Aacos[k L ■ y )} 2 . 

(2vr) Ji 


(2.14) 


Here D indicates that we have focused on contributions to the real-space power spectrum 
from observable wavemodes, and used Eq. (2.10) which is an expression valid for modes larger 
than /cl- The explicit y-dependence in the { } 2 parentheses represents the inhomogeneity of 
the power spectrum, including the hemispherical asymmetry. Note that the sinusoidal field 
fluctuation Eq. (2.9) sources power asymmetry not only with wavenumber /cl, but also with 
2 k L . 


3 Cosmological observables 

The best way to observationally constrain the power spectrum with broken translational 
invariance would be to do a Monte Carlo comparison with data using the spectrum in the 
form Eq. (2.10), which has not been attempted yet. Instead of doing that, we compute the 
multipole moments generated by the super-horizon fluctuations, in order to estimate the 
constraints on the power spectrum. 

Let us go back to the SAT expression (2.1) for ((x ), and substitute (2.9) for Aa(x). We 
represent the center of our observable Universe by xq and rewrite as 

k L ■ x = k L ■ (x - x Q ) + k L ■ x 0 , (3.1) 

also introducing 

6 = &l • xo , (3.2) 

as our particular phase in the large-scale mode, see Figure 1. For a specific observer 0 takes 
a fixed value, and our results will show how the observables depend on the location that the 
observer occupies on the super-horizon mode. As our observable patch, given by the distance 
to the last-scattering surface Xd ec , is taken to be much smaller than the wavelength of the 
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Figure 1. The observable Universe located in the super-horizon mode. The observable patch is 
centered around x 0 and has radius of order idee- We assume throughout \k^ ■ (x — x Q )\ -C 1 for cc 
within our observable patch. 


super-horizon mode, we expand Eq. (2.1) within our patch in terms of |■ (x — Xq)\ <C 1. 
This yields 

C( x ) =M a (dcj(a;) + Aa cos 9} + -M a u {5a(x) + A a cos 9} 2 


- {k L ■ (x - * 0 )} 


{Na + J\[ aa 8a(x)} A a sin 9 4— AT aa Aa 2 sin(20) 


; { k L • (x — So)} 2 [{Wcr + Afa<r8(j{x)} Aa COS 9 + Af aa A(J 2 cos(2#)] 

+ g {^L • (s - So)} 3 [{Wcr + A r ca 8a(x)} A a sin 9 + 2M aa Aa 2 sin(26*)] 
+ O {k L ■ (s - s 0 )} 4 + (terms with J\f arTrT , ■■■). 


(3.3) 


1 


Terms quadratic in k^ ■ (s — so) source quadrupole moments in the CMB, while cubic 
terms source octupole moments. The usual contributions to these from the continuum modes 
are given by the 5a terms in the first line, and can be computed once the spectrum Pg a is 
specified, while the contributions associated to the super-horizon mode are given by the 
other lines in this equation. The cross-term 5a(x)Aa in the first line represents the shift of 
the average field value within the observable patch from the background value of the entire 
universe, as we will soon see. The continuum modes contribute to each m-component of the 
multipoles, but by aligning our coordinate system to the direction of the super-horizon mode 
circular symmetry ensures that its contribution is only to the m = 0 mode of the multipoles, 
e.g. 020 for the quadrupole. 

A detailed derivation of the multipoles generated by super-horizon perturbations that 
can be expressed as a Taylor series, as in Eq. (3.3), was performed in Ref. [8]. Using their 
results, 3 the super-horizon mode contribution to the quadrupole and octupole are obtained, 
respectively, from the third and fourth lines of Eq. (3.3) as 

Aa2o ~ 0.10 x (k^Xdec) 2 {N a Aa cos 9 + N a oAa 2 cos(2#)} , (3-4) 

Aa 3 o ~ —0.019 X (k^Xdec) 3 {A/’ o -A<7sin0 + 2J\f CTa Aa 2 sin(26 1 )} . (3-5) 

3 In Ref. [8], multipole moments arising from a gravitational potential 'l' with a super-horizon mode were 

computed. Their results can be applied to the expression (3.3) for the curvature perturbation after the 
conversion T = — 3£/5. 
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Here we have dropped terms proportional to J\f a(7 da , supposing they are negligible compared 
to A f a terms. 1 Moreover, Xdec denotes the comoving distance to the surface of last scattering. 
The super-horizon mode contribution to the dipole cancels out at leading order in ^L^dec ^ e ’ 
tween the Sachs-Wolfe, integrated Sachs-Wolfe, and Doppler effects as was shown in Ref. [8] . 
Furthermore, since |aio| is not well constrained due to our peculiar velocity, we do not con¬ 
sider the dipole in the following discussions. We should also remark that the numerical factors 
in Eqs. (3.4) and (3.5) have dependence on the cosmological parameters. 

The terms quadratic in Act in Eqs. (3.4) and (3.5) arise from the non-linearities in 
the conversion process of field fluctuations into curvature perturbations. These effects were 
discussed in the context of the curvaton scenario in Ref. [8], and were also focused on in 
Ref. [14], 

The { } 2 parentheses in the real-space power spectrum (2.14) can also be expanded in 
terms of k^ ■ (y — yo), giving 


1 

(27r) 3 


'fc»fcL 


d 3 ke ik < x -^P ScT (k) {N a + N aa Ad cos 9f 


1-2 {k L -(y- y 0 )} 


A/’o-o-Acr sin 9 
N a + N a o Ad COS 9 


+ 0{k L ■ (y - yo)Y 


(3.7) 


Thus the all-sky averaged value of the power spectrum is obtained as 

P((k) = P$ a {k) (N a + A/^Ad cos 0) 2 

= (k) ■Mr|< rbg+ Aofcos0 ’ 


(3.8) 


where the vertical line indicates evaluation at the given field value. 4 5 By expanding this last 
expression in terms of A d, and dropping terms with N a aa, N a oaa, • • •, the second line reduces 
to the first. The final expression in Eq. (3.8) indicates that the average field value within our 
observable patch is shifted from the background value of the entire universe by Ad cos 9, at 
leading order in the &l expansion. We stress that this shift of the mean field value needs to 
be taken into account upon evaluating cosmological observables. 

The term linear in &l ■ {y — yo) in (3.7) gives the hemispherical power asymmetry. Half 
the difference of the power spectrum between spatial points that are separated by xaec^L/^L 
is 


m = 


1 A P ( {k) 

2 P((k) 


(/cLXdec) 

(/^L^dec) 


NooAd sin 9 
No- + NooAd cos 9 
NooAd sin# 


No 


crbg+Acr cos 6 


(3.9) 


whose observed value is 0.06 [6]. 

The asymmetry induced by the super-horizon mode is generated by the non-zero value 
of Noo ; which also generates cosmic non-Gaussianity. Hence the presence of the asymmetry is 


4 In the absence of the super-horizon field fluctuation Aa, then /nl 
yielding 

Moa5a(x) 6 . . 

AA "5 /nlC(x) - 


(5/6)(AA ct /A/^) and £(x) ~J\f a 6a(x), 

(3.6) 


When \((x)\ ~ 10 -5 , then as long as |/nl| <C 10 5 , the quantity (3.6) is tiny. 

’Here, strictly speaking the terms in the expansion (3.7) with even powers of (Icl • (y — yo)} further shift 
the averaged power spectrum, but we ignored such contributions as they are suppressed by powers of ^L^dec- 
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associated with non-Gaussianity and it is necessary to check whether the observed asymmetry 
is compatible with limits on non-Gaussianity. We introduce the non-linearity parameter in 
the observable patch, with phase 9, by 


/nl ( k ) 


5 A faa_ 

6 m 


fJbg+Acr cos 6 


5 _ Mgg _ 

6 + MggAa cos 9) 2 


(3.10) 


The middle expression in Eq. (3.10) is particularly interesting, since it lets us evaluate the non¬ 
linear parameter /nl in the asymmetric sky with the same expression that would be used in 
the absence of the asymmetry, but evaluated at a shifted value of the background field. That 
means that we can keep the usual expression for /nl, and evaluate it at the displaced value 
of the held, a = <7b g + Ad cos 9, according to our particular location in the long-wavelength 
perturbation. The same interpretation of a usual quantity evaluated at the displaced value 
of the background held also applies to the power spectrum expression in Eq. (3.8). It should 
be noted that the bispectrum in the presence of the super-horizon fluctuations should also 
exhibit broken translational invariance (as for the power spectrum (2.10)), and thus does not 
take conventional forms such as the local-type template. We use /nl defined in Eq. (3.10) 
simply as a rough estimate of the size of the bispectrum. 

Before ending this section, let us remark that we choose the initial hat hypersurface 
for the 5M calculations to be when a relevant wavemode exits the Hubble horizon during 
inflation. To be more specihc, we focus on a certain wavenumber k for which we consider 
the power asymmetry A(k), and define the time tk as when k = aH . The quantities in 
Eqs. (3.8), (3.9), and (3.10), such as P$g and Ad, are evaluated at t = tk, and M a , M aa , etc. 
are considered as derivatives in terms of a at t = tk- As for the multipoles, the quantities on 
the right-hand sides of Eqs. (3.4) and (3.5) can be evaluated at any time, because the initial 
hypersurface can be chosen arbitrarily for computing £ in Eq. (3.3). Hence, for convenience, 
we also evaluate the terms in Eqs. (3.4) and (3.5) at t = tk- 


4 Generic bound on the power asymmetry and /nl 


We now derive a bound on the power asymmetry in terms of the non-linearity parameter /nl- 
By combining the expressions for the quadrupole (3.4), power asymmetry (3.9), and non¬ 
linearity parameter (3.10), we obtain 

Aa 2 o /nl(&) ~ 0.09 x f3A 2 (k), where /3 =—(cos 29 + ~ cos fA . (4.1) 

sin 6 \ AlggAa J 


This shows that to generate a significant power asymmetry, one needs large enough values of 
both the (super-horizon mode induced) quadrupole and the non-Gaussianity, each of which 
is constrained by data. The value of j3 can be much smaller than unity only when the cosine 
terms are tuned to cancel each other — see Figure 2. For instance, when \M a /(MggAa)\ <C 1 , 
choices of 9 ss 7r/4, 37 t/ 4, • • • suppress (3 due to the quadrupole being suppressed, cf. Eq. (3.4). 
Otherwise, without fine-tuning of the model parameters, |/3| is of order unity or larger. Hence 
using |/31 > 1, we arrive at a generic bound 


|/nl(&)I >0.09x 


A 2 (k) 

|Aa2o| 


>30x 



(4.2) 


For the final expression we adopted a limit on the super-horizon mode contribution to the 
quadrupole of |Aa 2 o| 1.0 x 10~ 5 (c.f. Refs. [7, 11]); an exact constraint would require 
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e 

Figure 2. (3 as in Eq. (4.1) evaluated for Af a /(AT arT Acr) = 0.0033, 3, and 30 (blue, red, and green 
respectively). The dashed lines mark the regions where |/?| is or order unity. 


a model for the continuum contribution to the quadrupole but for our purposes requiring 
that this extra contribution is sub-dominant suffices. It is clear from Eq. (4.2) that, in 
order to realize A(k) ~ 0.06 as suggested by Planck [6], |/nl| has to be at least ~ 30. The 
required /nl for producing the power asymmetry becomes even larger when \(3\ 3> 1; then, to 
be compatible with the current constraints on local-type non-Gaussianity, /nl = 2.7T5.8 [22], 
the non-linearity parameter needs to possess a scale-dependence. 

5 Example parameters 

In this section, we present some example parameter sets that allow hemispherical power 
asymmetry at the level suggested by observations without violating other constraints. Such 
a large asymmetry is generically accompanied by the non-linearity parameter |/nl| at low 
is being larger than ~ 30, cf. Eq. (4.2). Unless this lower bound on /nl is circumvented by 
fine-tuning of the phase 9 to suppress f3, a running of /nl is needed to accommodate the 
CMB constraints on non-Gaussian perturbations, but anyway the asymmetry is observed 
to be scale-dependent, dying away at high multipoles [2, 4, 23, 24], and hence one should 
anticipate running of the non-Gaussianity in a viable model. 

Let us start by summarizing the constraints that should be satisfied. Firstly, we have 
been assuming that the power asymmetry is seeded by a field fluctuation with wavelength 
much larger than our observable Universe, 

khXdec < 1- (5.1) 

Furthermore, upon using the 5J\f formalism, we have considered the Universe to be homoge¬ 
neous and isotropic on very large scales, thus 


\Af a Aa\ , |A/' (7(T Acr 2 


- 9 - 


< 1. 


(5.2) 







It should also be noted that we have assumed 






-Sa(x) 


< 1 


(5.3) 


when computing the multipoles (3.4) and (3.5). The averaged power spectrum 

2ir 2 

Pdk) = (5-4) 

whose expression is given in (3.8), needs to satisfy the Planck normalization [25] 

V c (k) « 2.2 x 10" 9 . (5.5) 

For the quadrupole (3.4) and octupole (3.5) induced by the super-horizon field fluctuations, 
as described above we adopt the constraints 

|Aa 2 o|,|Aa 3 o| < 1-0 x 10- 5 . (5.6) 


In principle, the induced multipoles may partially cancel out the multipoles from the con¬ 
tinuum spectrum (i.e. perturbations in the absence of super-horizon fluctuations), and then 
Eq. (5.6) would not apply, but this would be less likely than the original multipoles being 
coincidentally small. We do not consider such cases in the following discussions. 

The observables such as V^(k), A(k), and Aa 2 o are expressed in terms of five quantities 

Af a Aa, J\f aa Aa 2 , k L x dec , , 0, (5.7) 

where Vs a is related to the power spectrum of the field fluctuations P§ a as in Eq. (5.4). In the 
following, we look into two cases, where Af a Aa is larger or smaller than A^o-Acr 2 respectively. 
Let us repeat that, unless indicated otherwise, M a and M aa are global universe background 
values, i.e. quantities evaluated at 


5.1 Cases with |A/’ (J Acr| ~ \N aa A(j 2 1 and |A/’ (J Acj| 2> \N au A(j 2 1 
In this case, the amplitude of the power asymmetry is given by 

A{k) ~ -(fcLffdec) - ^ 7 ^- sin 9, (5.8) 

unless A/’o-Ad and A/’o-o-Ad 2 happen to cancel each other in the denominator of Eq. (3.9). 
From the assumption \J\T a Aa\ > |A/’ (7(7 Acj 2 | , one sees that the wavelength of the super-horizon 
fluctuation cannot be much larger than the size of the observable Universe if it is to seed a 
power asymmetry as large as |A(/c)| ~ 0.06. 

The lower bound on /nLj Eq. (4.2), is saturated when the parameter f3 in Eq. (4.1) is 
of order unity. In particular, when |A/’ 0 -A(t| |A/’ cr(J A<7 2 |, our observable Universe needs to 

be located at a special position along the super-horizon sinusoidal fluctuation such that cos 9 
nearly vanishes to have |/3| < 1, cf. Figure 2. 

An example set of parameters that gives the observed |A| while |/nl| 10, enabled by 

tuning the phase 9 , is 


A/’o-Acr = 0.030, N a( jA(j 2 = 0.010, /q^dec = 0.20, 


Vsa(k) 

A a 2 


= 2.1 x 10 ~ 6 , 


9 = 1.3, 


(5.9) 
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leading to a suppressed /? ~ —0.06, and yields 


V c (k) « 2.2 x 10- 9 , A(fc)«-0.06, /nl(&) ~ 8, 

c « (5.10) 

Aa 20 ~ -2.3 x 10' 6 , Aa 30 »-6x 1(T 6 . 

Furthermore, if the typical amplitude of 6a(x) is given by the field’s power spectrum as 
1 /2 

5a(x) ~ V s (k), then the constraint Eq. (5.3) is also satisfied. 

5.2 Case with \M a Aa\ <C |A/’ CTC rA(j 2 | 

Here we study the limit considered by Kanno et al [14], which focused on the non-linearity 
in the field fluctuations. When the quadratic contribution A/" 0 - cr Ad 2 dominates over the 
linear A/" ct Act, then unless cos0 is close to zero, the power asymmetry is directly tied to the 
wavenumber of the super-horizon mode, 


A{k) ~ -(fc L x d ec)tan6», 


(5.11) 


and (5 is set merely by the phase, i.e., f5 ~ cos 20/ sin 2 6. In this limit, our relation (4.1) 
simplifies as (ignoring phase factors, hence use of ~), 


|Aa 2 o /nl| ~ 3 x 10 4 



(5.12) 


agreeing with Eq. (7) of Ref. [14]. 6 A similar expression was also obtained in the context of 
the curvaton scenario in Eq. (10) of Ref. [7], though it should be noted that the parameter A 
in Ref. [7] is defined as A P^/P^, i.e. twice the A used in this article. 

Now, \/3\ < 1 is realized in a rather large range of phase values around 0 ~ 7t/2 and 
37t/2, as can be seen in Figure 2. For instance, a set of parameters that saturates the lower 
bound (4.2) is 


giving (5 


AT^Aa = 0 . 0010 , 

VsAk) 

Act 2 


A/’ (7( jAct 2 = 0.15, k^x^ec = 0.040, 
= 3.3 x 10 -7 , 0 = 1.0, 


—0.6, and 


Vc_{k) « 2.2 x 10~ 9 , A(k) « -0.06, / NL (A:) « 19, 

Aa 2 o ~ —1.0 x 10 Aa 3 o ~ —3 x 10 ". 


(5.13) 


(5.14) 


1 /2 

If 5(j{x) ~ V Sa (k) further holds, then the constraint Eq. (5.3) is also satisfied. 

However, the condition |A/’ (J Act| <C |A/’ (J(T Act 2 | may be suggesting the breakdown of the 
Taylor expansion of A f for a = dbg ± Act. If this is actually the case, then the computation 
ignoring terms with M a aa and higher-order derivatives is invalid, motivating a more general 
calculation that we undertake in the next section. 

6 Ref. [14] studied a case where two fields contribute to the curvature perturbations. However, there it was 

assumed that only one of the two fields gives non-linear contributions, and thus their Eq. (7) matches with 
our single-field result. 
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6 Full account of all non-linearities 


Until now we have expanded (derivatives of) the e-folding number J\f in terms of Acr(x) to 
carry out the integrations. However the expansion can break down when the super-horizon 
field fluctuation Act is much larger than the small-scale fluctuations 5a, and there does not 
exist a hierarchy of \Af a Aa\ |A/’ (7(T A3 :2 | |A/’ cr(TO -Ac7 3 ^> ■ ■ ■. Furthermore, we saw in Sec¬ 
tion 5 that such a hierarchy is typically absent in parameter regimes capable of explaining the 
observed asymmetry. In this section we generalize our treatment by avoiding Act expansions, 
instead expanding quantities directly in terms of the wavemode of the super-horizon fluctu¬ 
ations &l (in units of the inverse of the horizon distance cc^ec)- The k\, expansions should be 
valid whenever considering fluctuations with wavelengths much larger than the size of the 
observable Universe. 

The curvature perturbation is expressed in terms of the 5.\f formalism as 

C(*) = a(x) - ^\a bs > (6T) 

where the field value at each spatial position on the initial flat hypersurface is written as 


a(x) = <7bg + 5a(x) + Aa(x). 


( 6 . 2 ) 


While Act may be large, we consider the perturbative expansion in terms of the continuum 
spectrum of fluctuations 5a to be valid. Thus we first expand ( in terms of 5a, 


C(x) - J\f\ CTbg+ ACT(a:) ^\cT b g 

+ 5a(x) A 


+ 6(j2 W ^( 2 ) 

cr bg +Acr(x) 2 


a b „+Au(x) 


+ 


(6.3) 


where the nth derivative of A f in terms of the field a is denoted as J\f in this section. Note 
that the e-folding derivatives are evaluated at a field value shifted by Act. 

As in the earlier treatment, to compute the power spectrum we treat the fluctuation Act 
as a non-stochastic quantity with the sinusoidal form, 


A a(x) = Act cos (k l • x ), (6-4) 

and also adopt the correlation relations (2.4), (2.5), and (2.6) for 5a. Then the power 
spectrum of the Fourier modes of £ can be written as 


(C fcl Cfc 2 ) = J d 3 *d 3 2/e- i(fei a:+fc2 y) 






a bs +Air(x) ~ yV U bg ; \a bg +Aa{y) ■ • Kr bg 

+ 2 W ( } ffbg+A(T(y) ( S(7 (l/)> 


+ 2 (^Ug+A a(y) N \* bg ) ^ ] abs+Aa(x) (&<* (*)> 


(6.5) 




(5a(x) 5a(y)} + 0(5a) 4 


o- bg +A(j(3:) 

Obg+A cr(y) 
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Now, unlike in the previous sections, we directly expand quantities in terms of k^ ■ (x — 
xq), without carrying out Act expansions. This gives, for the derivatives of Af, 


A7 (n) 


rr bg +Acr(x) 


= Afn'^ — &l • (* — *o) ATn n+1 ' > Act sin 6 


{fcL ' ( x ^o)} / i r(n+l) 


{n+2) Aa 2 S m 2 9 


(n+3) 


Act 3 sin 3 6 


Afg Ad cos 0 — M g 

+ — ——^— ir °^ (^J\fg n+1 ' > Act sin 6 + 3Afg U+2 ^ Ad 2 sin 0 cos 9 — Af' t 
+ 0{k L • (s-cco)} 4 . 

Here Afg is an abbreviation for 

Af^ n) = M {n) 

(Tb g +A ct cos 6 

which is an ^-independent quantity. Expression (6.6) can be Fourier expanded using 
1 


( 6 . 6 ) 


d A x e lkx {&L • (x — x 0 )} n = e 


(27r) s 


n _ —ik-XQ 


(' ik L -V k ) n d(k ), 


(6.7) 


( 6 . 8 ) 


yielding 


^(n) _ 1 

yVfc “ (2 tt)3 


d 3 x e~ ik ' x M (n) 


(T bg +A a(x) 

= Af^ n) S(k)-Af^ n+1) Adsmde- ik X0 ik L • V fc 5(fe) 


(6.9) 


" +1 ^ Act cos 0 — J\rf} n+2 ' > Act 2 sin 2 o) e lk ' x °{iki J ■ Vfc) 2 5(fc) + 0{ki) 


Here we have denoted terms that contain cubic or higher orders of k l by 0 (&l) 3 • In terms 
of A/l ?t \ the power spectrum (6.5) can be written as 


(CfciCfca) = (2^) 6 {a4 x - AA| ffbg h(fc!)} {a4 2 - 5(k 2 )} 


+ 


(2vr) £ 


rf 3 P^(ri [a/£ } {A/- fe2 - A7| CTbg <S(fc 2 )} +xg {a/ - * - A/| CTbg <$(fci)}] (6.10) 


+ (2vr) 3 y d 3 p^(p)AAW_ p <l P + °(^) 4 - 

Here we note that Eq. (6.8) vanishes for k ^ 0, and so does Aup to finite order in k l- 
Thus, focusing on nonzero wavenumbers k±,k 2 ^ 0, the first and second lines on the right- 
hand side of Eq. (6.10) can be dropped. After carrying out integration by parts on the third 
line, and ignoring the scale-dependence of the power spectrum of the held fluctuations by 
assuming 

0 = (ik L ■ V p ) m P Sa (\k - p\)\ Q for m = l,2, •••, (6.11) 
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one finds 


(CfciCfc2 


(2vr) 


{P$ a (kl) + PsAfa)} 


K 


( 1 ) 


m 


- 2A/"i 1) A/ r i 2) A5 sin 9 e~ ik x ° (zfc L ■ V fc )<5(fe) 


— i Ng 1 '*Nq 2 ' 1 Act cos 


«-(<X <S, + (Aff) 2 ) 


Act 2 sin 2 0 


( 6 . 12 ) 


xe-‘^(*fe L -V fc ) 2 i(fc) 


k=ki+k2 


+ 0(k L ) 3 + 0(5a) 4 . 


The terms with derivatives of the delta function represent the broken translational invariance. 

Hereafter we neglect terms of O(da) 4 and 0(k\ J ) 3 . Also ignoring contributions from 
zero wavemodes, the power spectrum in real space is computed by integrating Eq. (6.12), 
giving 


(a x )C ( y ))=— 3 J d 3 ke ik < x -^PsAk) 


A/g 1 ^ — 2k l • (y — ®o) JVf} 1 ' 1 JVg 2 ' 1 Aa sin 9 


- {k L ■ ( y - xq)} 2 |a/' 0 (1) A/' 0 (2) Actcos(9 - + (a/^) ^ 


Act 2 sin 2 6 


(6.13) 


One can read off the all-sky averaged power spectrum, at the leading order in &l, 

Pdk) = Psa(k) (a/^) 2 , 

and the hemispherical power asymmetry 

_1A P c (k) _ AT^AasmO 

2 P c (fc) lfcLXdecj jtfCD • 

u 

We also define the non-linearity parameter as 



/nl(&) 


5 A/f 

6 K ) ) 2 ' 


(6.14) 


(6.15) 


(6.16) 


The multipoles that arise from the super-horizon fluctuations can be computed by expanding 
C(x) in terms of k^-(x-xo) (cf. Eqs. (6.3) and (6.6), and see also discussions around Eqs. (3.4) 
and (3.5)), 


Aa 2 o ~ 0.10 X ( ki,Xdec ) 2 (A/’g 1 ) Act cos 9 — Afjj^Aa 2 sin 2 9^ , (6-17) 

Aa 3 o ~ —0.019 x (/cL^dec) 3 Ad sin 9 + 3J\Tq^ Act 2 sin 0 cos 0 — A/)) 3,1 Act 3 sin 3 9^ . (6.18) 
Combining the above expressions gives 

Aa 20 / N L(A;)«0.09x A 2 (fc) (-1 +^. (6.19) 

V Nq Ad s m 2 9 J 
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Further introducing 


5Nl(&) 


25 

54 



3 


to parameterize the size of the trispectrum, one also obtains 


Atl30 


/nlW 

9NL{k) 


0.023 x A 3 (k ) -1 + 


3Acos 9 

A 3 ) 


+ 


K 


(i) 


Aff’ Act sin 2 9 Afi 3 ' ) Act 2 


sin 


( 6 . 20 ) 


( 6 . 21 ) 


though its implication is less clear compared to Eq. (6.19). 

The expressions in the previous sections are reproduced by expanding the formulae in 
this section by A a and ignoring terms of A, AA 4 ), • • •. One can check that the expressions 
(6.13), (6.14), (6.15), (6.16), (6.17), (6.18), (6.19) reduce to, respectively, (3.7), (3.8), (3.9), 
(3.10), (3.4), (3.5), (4.1). 

Even though we have not carried out Act expansions in this section, most of the results 
contain terms only up to cubic order in A a. This is because Ad is generically accompanied 
by the wavenumber &l, and thus terms with higher orders of Ad are suppressed by higher 
powers of k^. 


7 Summary 

We have presented a comprehensive analysis of curvature perturbations induced by field 
fluctuations with wavelength larger than our observable Universe. Treating the super-horizon 
fluctuations as non-stochastic, we evaluated correlation functions of perturbations on small 
scales. In other words, we have analyzed conditional probability distributions of small-scale 
perturbations given a fixed form for the super-horizon fluctuation. 

Our formalism improves on previous approaches by consistently incorporating the super¬ 
horizon and small-scale fluctuations using the 5Af formalism. In particular, we obtained 
the full expression for the power spectrum of the curvature perturbations that manifests 
translational-invariance breaking in the observable patch, Eqs. (2.10) and (6.12). We also 
point out the importance of choosing the appropriate ‘background’ value where various quan¬ 
tities should be evaluated, which leads to significant simplifications in the final expressions. 
Observational signatures that accompany the power spectrum with broken translational in¬ 
variance were computed, such as the hemispherical power asymmetry and multipole moments 
in the CMB, including the full dependence on the phase 9 of our location within the super¬ 
horizon mode. 

Non-linear effects are introduced through the conversion process of the field fluctua¬ 
tions into curvature perturbations, and also as the perturbations source the CMB temper¬ 
ature fluctuations. As a consequence, the observational signatures show up not only at the 
wavemode of the super-horizon field fluctuation, but also at smaller length scales. In the 
first half of the article, we carried out analyses by expanding various quantities in terms of 
the wavenumber k l and amplitude Acf of the super-horizon field fluctuation, and dropped 
terms with cubic or higher orders in Air. We then further generalized the calculations in 
Section 6 by carrying out expansions only in terms of £;l; but not in terms of Acf, in order 
to accommodate potentially large higher-order derivatives, enabling calculations such as the 
trispectrum of the curvature perturbations. 
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Our formalism can be applied to any model where the curvature perturbations are 
sourced by a held that has super-horizon fluctuations. We have shown that there exist 
parameter ranges where the hemispherical power asymmetry can be as large as the observed 
value A ~ 0.06, while the non-Gaussianity is suppressed to |/nl| 1$ 10- It is important to 
construct explicit models that realize such parameter ranges. Alternatively, super-horizon 
fluctuations in models with running /nl may also realize significant power asymmetry and still 
be consistent with existing constraints on non-Gaussianity. We will report on investigations 
of explicit models in a future publication. 

We have obtained constraints on primordial power spectra with broken translational 
invariance mainly through discussing the induced multipoles in the CMB. However, a direct 
way to constrain such power spectra would be to do a MCMC comparison with observational 
data using the full expression of the spectrum presented in Eq. (2.10). It would be interesting 
to see how the constraint on the power asymmetry changes, compared to current limits 
imposed without having a model for the power spectrum. The calculations in this article 
can be extended to compute higher-order correlation functions as well. The higher-order 
functions should also exhibit translational invariance breaking, giving rise to new shapes for 
the bispectrum and trispectrum. 

Spatial dependences of statistical quantities can, albeit indirectly, give us hints about 
cosmological structures beyond our observable Universe. Our formalism is useful for probing 
super-horizon perturbations from observations of asymmetric signals in the sky. 
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